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Lorentz gauge fixing and the Gribov problem: 

the fermion correlator in lattice compact QED with Wilson fermions * 
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For the Lorentz gauge the influence of Gribov copies on the fermion propagator is investigated in quenched 
lattice compact QED. In the Coulomb phase zero-momentum modes of the gauge fields are shown to be the main 
reason for a significant deviation from ordinary perturbation theory. 
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Numerical lattice computations of gauge- 
dependent gauge or fermion field correlators can 
give us detailed information about nonperturba- 
tive properties of quantum fields and allow a di- 
rect comparison with perturbation theory in the 
continuum. However, most of the known iterative 
gauge fixing procedures are faced with the exis- 
tence of Gribov copies. For compact lattice QED 
within the Coulomb phase the standard Lorentz 
gauge fixing causes the transverse photon corre- 
lator to deviate significantly from the expected 
zero-mass behaviour [|l|,|[. Also the fermion cor- 
relator depends strongly on the achieved gauge 
copies 1^. Numerical [@,§-^ and analytical 
studies have shown that the main gauge field ex- 
citations responsible for the occurence of disturb- 
ing gauge copies are double Dirac sheets (DDS) 
and zero-momentum m,odes (ZMM). 

In the present talk we are going to discuss the 
influence of the ZMM on the fermion correlator. 
We shall demonstrate that only a proper account 
of them will allow us to determine the renormal- 
ized fermion mass. 

We consider 4d compact QED in the quenched 
approximation on a finite lattice {V — N^ x 
Nt). We employ the standard Wilson (plaquette) 
gauge and fermionic action, respectively. The lat- 
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ter is given by 

Sf = Y,i^Mxv{0)4,y, m = i-KD, (1) 
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where P^ = i±7M- The Ux,i_i — e'^^ ^' denotes the 
gauge link degrees of freedom. Periodic boundary 
conditions (b.c.) are implied. The fermion fields 
are anti-periodic in 2:4. The fermion correlator 
for a given gauge field 9 reads 

r(-'^) = ^EEM.-,U,,+.W. (2) 
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We shall restrict ourselves to the vectorial part 
1_ 
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-ReTr(74r(r;0)), 



(3) 



with the trace taken with respect to the spinor 
indices. 

Lateron, we shall compare the average ( Fy )e 
within the approximation where only constant 
gauge fields are taken into account. The corre- 
lator in a uniform background 9x.p, = (t)^,, — vr < 
0^1 < TT, M = 1. ■ • • : 4 is given by 
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2 (l-t- £2JVt _ 2£'^* cos{(l>4Nt)) 
X ([£■"+ £2A'*-]cos(04t) (4) 
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For (p^ — 0, fi — 1,- ■ ■ ,4 the free ferniion 
correlator for finite lattice size g is reproduced. 
Note, that the bare mass toq is related to the 
hopping-parameter by k = 1/(8 + 2too). 

In numerical simulations the Lorentz gauge is 
fixed by iteratively maximizing the functional 
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with respect to (periodic) gauge transformations. 

The algorithm is called standard Lorentz gauge 
fixing, if it consists only of local maximization 
and overrelaxation steps. It is well-known that 
this procedure gets stuck into local maxima of the 
functional (S) (gauge copies). It is a common be- 
lief that the Gribov problem has to be solved by 
searching for the global maximum providing the 
best gauge copy. In we have shown that in 
order to reach the global maximum we have nec- 
essarily to remove both the DDS and the ZMM 
from the gauge fields. 

DDS can be identified as follows. The plaque- 
tte angle dx,)!^ is decomposed into the gauge in- 
variant (electro-) magnetic flux dx.tiu S ("""i ""] 
and the discrete gauge-dependent contribution 
27rna;^^y, Ux^fjiu = 0, ±1,±2. The latter repre- 
sents a Dirac string passing through the given 
plaquette if nx,fj,u = ±1 {Dirac plaquette). A 
set of Dirac plaqucttes providing a world sheet 
of a Dirac string on the dual lattice is called 
Dirac sheet. DDS consist of two sheets with op- 
posite flux orientation extending over the whole 
lattice and closing themselves by the periodic b.c. 
They can easily be identified by counting the total 
number of Dirac plaquettes Nf^p for each choice 
{lJ,\v)- DDS can be removed by periodic gauge 
transformations. For standard Lorentz gauge fix- 
ing DDS have been seen to occur independently 



of the lattice size and /3 yj5|,|6|] . As a consequence 
the non-zero momentum transverse photon cor- 
relator significantly deviates from the expected 
zero-mass behaviour [gj-Q]. 
The ZMM of the gauge field 
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do not contribute to the pure gauge field action ei- 
ther. For gauge configurations representing small 
fluctuations around constant modes it is easy to 
see, that the global maximum of the functional 
(|5|) requires 0^ = . The latter condition can be 
achieved by non-periodic gauge transformations 
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We realize a proper gauge fixing procedure as 
proposed in Q . The successive Lorentz gauge it- 
eration steps are always followed by non-periodic 
gauge transformations suppressing the ZMM. Ad- 
ditionally we check, whether the gauge fields con- 
tain yet DDS. The latter can be excluded by re- 
peating the procedure with initial random gauges. 
We call the combined procedure zero-momentum 
Lorentz gauge (ZML gauge). It provides with 
very high accuracy the global maximum of the 
gauge functional. The photon propagator per- 
fectly agrees with the expected perturbative re- 
sult throughout the Coulomb phase M. 

Our Monte Carlo simulations were carried out 
with a filter heat bath method. In order to extract 
the pure ZMM effect, we first apply the standard 
Lorentz gauge procedure modified by initial ran- 
dom gauges in order to suppress DDS. We de- 
note this modified standard Lorentz gauge proce- 
dure by LG. We compare the latter with the ZML 
gauge described above. For both these gauges 
we have computed the averaged fermion correla- 
tor employing the conjugate gradient method and 
point-like sources. In the upper part of Fig. yj we 
have plotted Tv{t;0) (normalized to unity at 
T = 1). The situation seen is typical for a wide 
range of parameter values within the Coulomb 
phase. Obviously, there is a strong dependence of 
the fermion propagator on the gauge fixing proce- 
dure resulting in the presence or absence of ZMM. 
The masses to be extracted seem to have differ- 
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Figure 1. The fermion vector propagator (a) and 
the effective mass (b) at /3 = 2 and k = .122 
on a 12 X 6'^ lattice for LG and ZML gauges as 
explained in the text. 



Figure 2. Free fermion propagator (dashed hne) 
and averaged constant-mode propagator in the 
LG case (full line) for /? = 2, 10, k = .120, .124, 
lattice size 12 x 6^, 16 x 8^. 



ent values. Let us determine the effective mass 
meffir) in accordance with 



(r(r + l;g))e _ cosh[^(T)(JVt/2 - 
(r(r;e))e 



1)] 



cosh[£;(T)(iVt/2-T)] 



(7) 



where E{t) = ln(TOe//(T) + l). See the lower part 
of Fig. 1^. In the LG case no plateau is visible, 
whereas the ZML case provides a very stable one. 
Thus, only the ZML gauge yields a reliable mass 
estimate, whereas the standard method to fix the 
Lorentz gauge obviously fails. To get deeper in- 
sight into the effect of ZMM for the LG case (with 
DDS suppressed) we measure the probability dis- 
tributions P{4>) for the space- and time-like com- 
ponents of ZMM according to Eq. (^. The dis- 
tributions turn out to be flat up to an effective 
cutoff at |0;j| — T^/Nfj. and to be independent of 
p. In accordance with Eq. (^ we compute the 
fermion propagator for constant modes in the LG 
case and average 



(ry(T;0))^= /[d0]P(0)ry(T;0). 



(8) 



The results for several parameter sets are pre- 
sented in Fig. together with the corresponding 
free (i.e. zero-background) propagator. 

The model describes qualitatively the influence 
of ZMM on the (full) fermion propagator in the 
LG case very well. It shows that the ZMM effect 
should not be expected to become weaker with 
increasing /3 and/or lattice size. 
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